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ABSTRACT 


A  design  for  a  moving  target  indication  (MTI)  subsystem  for  the  infrared 
airborne  radar  (IRAR)  is  investigated.  A  Doppler  MTI  system  incorporating 
heterodyne  detection  and  a  cw  C02  laser  is  determined  to  be  the  most  suitable 
approach.  Numerous  analyses  of  topics  including  field-of-view  coverage, 
Dopper  resolution,  signal -to-noise  ratio,  clutter  effects,  etc.  are  presented 
which  serve  to  define  component  specifications  and  processing  requirements. 
The  potential  architecture  of  the  MTI  system  is  discussed  and  the  impact 
of  developing  technologies  on  system  and  component  design  is  estimated. 
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1.  INTRODUCTION 

1-3 

The  goal  of  the  infrared  airborne  radar  (IRAR)  program  is  to  develop 
a  fire  control  system  for  tactical  aircraft  utilized  in  a  close  air  support 
role.  The  concept  of  the  IRAR  system  calls  for  a  laser  moving  target  indica¬ 
tion  (MTI)  system  to  detect  targets  of  interest  (i.e.,  moving  targets)  and 
provide  a  location  cue  to  a  high-resolution  imaging  infrared  radar.  The 
imager  would  be  used  first  to  positively  verify  the  identity  of  a  suspected 
target  and  later  to  act  as  a  gunsight  for  the  aircraft's  ordnance  system. 

In  this  report  we  will  examine  the  design  criteria  for  the  IRAR  laser  MTI 
subsystem  and  identify  areas  which  require  further  investigation 

MTI  systems  can  be  divided  into  two  general  classes:  delay  line  cancel¬ 
lation  techniques  and  Doppler  techniques.  In  the  first  class  of  techniques 
a  radar  return  from  the  scene  obtained  with  one  pulse  is  subtracted  from  the 
return  obtained  from  the  preceding  pulse.  As  a  result  stationary  targets  will 
be  cancelled  out  and  moving  targets  will  be  emphasized.  This  form  of  MTI 
is  not  particularly  suitable  for  our  use  for  two  reasons.  First,  the  air¬ 
craft's  own  motion  will  cause  every  object  in  the  scene  to  shift  its  position 
making  it  difficult  to  detect  those  objects  which  are  truly  moving.  Second, 
delay  line  techniques  require  multiple  observations  of  the  scene.  Tactical 
aircraft  applications  may  not  permit  more  than  one  observation  of  each  point 
in  the  scene  as  the  aircraft  flies  along.  Consequently,  the  IRAR  system 
must  utilize  a  Doppler  MTI  technique.  In  a  Doppler  technique  the  frequency 
shifts  between  the  transmitted  and  received  signals  are  analyzed  to  deter¬ 
mine  the  relative  velocities  between  objects  in  the  scene  and  the  radar.  In 
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a  fixed-base  installation  any  Doppler -shifted  return  indicates  a  moving  tar¬ 
get.  In  an  airborne  application  any  difference  between  the  Doppler  shift  due 
to  the  relative  speed  of  the  aircraft  with  respect  to  the  ground  and  the 
Doppler  shift  of  an  object  on  the  ground  indicates  that  the  object  is  moving. 

Several  microwave  radar  MTI  systems  are  currently  in  use  on  tactical 
aircraft  and  could  in  principle  be  employed  for  cuing  the  IRAR  imager.  We 
have  opted  for  a  laser  radar  MTI  system  because  it  offers  a  possibility  for 
sharing  optical  components  with  the  imager  (and  thus  achieving  additional 
compactness)  which  a  microwave  system  does  not  have. 

The  basic  concept  of  the  proposed  IRAR  MTI  subsystem  is  as  follows .  A 
laser  beam  illuminates  a  patch  on  the  ground  ahead  of  the  aircraft.  A  com¬ 
bination  of  the  aircraft's  forward  motion  and  a  rapid  side-to-side  scanning 
of  the  illuminated  patch  serves  to  illuminate  the  entire  ground  (at  least 
once)  within  a  specified  search  field-of-view.  The  laser  radiation  reflected 
by  the  scene  is  collected  by  a  telescope  and  imaged  onto  a  detector  system 
(possibly  consisting  of  an  array  of  detectors) .  The  received  radiation  is 
frequency  analyzed  to  yield  the  relative  radial  velocity  between  the 
aircraft  and  each  portion  of  the  scene.  Comparison  with  the  aircraft's 
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ground  velocity  (possibly  obtained  by  the  IRAR  MTI  system  itself  or  possibly 
provided  by  another  aircraft  subsystem)  yields  those  portions  of  the  scene 
which  are  moving  with  respect  to  the  ground.  The  MTI  system  then  directs 
the  IRAR  imager  to  image  the  moving  areas  for  possible  hostile  target  identi¬ 
fication. 

Although  it  is  highly  desirable  to  have  the  IRAR  MTI  subsystem  share  as 
many  components  as  possible  with  the  IRAR  imager  subsystem,  we  will  not 
restrict  ourselves  by  arbitrarily  imposing  this  restraint.  We  begin  the 
design  analysis  by  corparing  the  relative  merits  of  various  detector  systems 
and  laser  systems.  These  comparisons  dictate  the  use  of  a  cw  002  laser- 
heterodyne  detector  system  to  perform  the  MTI  function.  Next  the  conditions 
imposed  by  the  need  for  efficient  heterodyne  detection  and  adequate  coverage 
of  the  search  field-of-view  on  the  size,  shape,  and  structure  of  the  hetero¬ 
dyne  detector  system  are  investigated.  These  aspects  of  the  problem  dictate 
that  a  large  one-dimensional  array  must  be  employed.  Given  the  general 
design  of  an  MTI  system  as  determined  from  the  analyses  in  the  preceding 
sections,  the  signal-to-noise  ratio  (SNR)  for  realistic  values  of  the  param¬ 
eters  envisioned  in  actual  use  situations  is  estimated  and  the  beneficial 
effects  of  multiple  resolution  element  integration  are  investigated.  This 
analysis  indicates  that  resolution  element  integration  is  required  to  extend 
the  operational  regime  of  the  IRAR  system  to  include  moderately  bad  weather. 
The  effects  of  clutter  are  also  considered  and  clutter  is  found  to  have  a 
significant  effect  on  the  minimum  detectable  target  velocity.  Next,  the 
architecture  of  the  Doppler  frequency  analyzer  and  resolution  element 
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integration  systems  are  discussed.  The  potential  utility  of  various  emerging 
technologies  (acoustooptic,  surface  acoustic  wave  (SAW),  and  charge -coupled 
device  (CCD))  is  also  considered.  Throughout  this  report  the  compatibility 
of  the  total  MU  subsystem  resulting  from  this  work  with  the  proposed  IRAR 
imaging  subsystem  is  analyzed.  Although  the  electronic  processing  portions 
of  each  subsystem  are  obviously  different,  it  is  found  that  a  single  optical 
system  (laser,  transmit  and  receive  optics,  and  detector)  can  be  made  to 
suffice  for  both  the  imaging  and  MTI  functions. 

2.  CHOICE  OF  LASER  AND  DETECTOR  SYSTEMS 

In  the  Introduction  we  gave  arguments  supporting  the  use  of  a  Doppler 
MTI  technique  in  the  IRAR  system.  There  are  two  conceptually  different 
approaches  to  Doppler  MTI.  In  the  first  approach  the  laser  radiation  return¬ 
ing  from  an  element  of  the  scene  is  frequency -analyzed  by  a  spectrometer,  a 
set  of  narrow-band  filters,  or  some  other  optical  device.  The  dispersed 
optical  signal  is  then  converted  to  a  set  of  electrical  signals  by  direct 
detectors,  a  different  detector  being  used  for  each  spectral  resolution 
element  of  interest.  The  electrical  signals  are  subsequently  analyzed  to 
determine  the  Doppler  shifts  in  the  scene  from  which  the  existence  or  non¬ 
existence  of  moving  targets  is  deduced.  In  the  second  approach  the  returning 
laser  radiation  and  radiation  from  a  local  oscillator  laser  are  simultaneously 
detected  by  a  single  heterodyne  detector.  The  electrical  output  from  the 
heterodyne  detector  contains  frequency  components  which  are  directly  related 
to  the  Doppler  shifts  in  the  scene.  By  passing  the  signal  through  a  bank  of 
parallel  electrical  bandpass  filters  or  through  a  real-time  Fourier  transformer 
a  nunber  of  electrical  signals  are  obtained  which  can  be  processed  as  in  the 
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first  approach  to  determine  the  Doppler  shifts  and  thereby  deduce  the  exis¬ 
tence  or  non-existence  of  moving  targets. 


The  IRAR  MTI  subsystem  should  be  capable  of  2  kmph  motion  discrimination. 

To  accomplish  this  the  Doppler  MTI  system  must  have  a  frequency  resolution 

O 

capability  R'  =  v/6v  of  at  least  5.4  x  10  .  Here  <$v  is  the  minimum  resolvable 
Doppler  shift  imposed  on  a  transmitted  laser  beam  of  frequency  \>.  Almost  all 
optical  frequency-analyzing  devices  have  resolutions  which  fall  well  short  of 

4 

this  requirement .  For  example,  narrow-bandpass  optical  filters  have  R'  <  10  ; 
typical  grating  spectrographs  have  R'  ^  105;  and  a  very  good  Fabry-Perot  inter- 

g 

ferometer  may  have  a  value  for  R'  approaching  10  in  the  visible  region. 

Although  one  could  in  principle  build  an  optical  system  with  the  required 
resolution,  it  would  probably  be  quite  large,  very  conplicated,  and  extremely 
sensitive  to  environmental  factors  (vibration,  temperature,  etc).  Consequently, 
we  must  conclude  that  the  first  approach  is  unsuitable  for  use  in  the  IRAR 
system.  The  maximum  theoretical  resolution  of  the  heterodyne  detection 
approach  is  limited  only  by  the  laser  and  local  oscillator  frequency  stabili¬ 
ties  and  the  bandwidth  of  the  electrical  circuitry.  The  electrical  bandwidths 
can  be  made  as  narrow  as  one  desires  and  lasers  have  been  built  with  frequency 

stabilities  v/6v  which  are  orders  of  magnitude  in  excess  of  the  required  value 
0 

of  5.4  x  10  .  Thus  there  is  no  fundamental  limitation  to  using  the  heterodyne 
approach  in  the  IRAR  MTI  subsystem. 

The  signal -to -noise  ratio  obtained  from  a  heterodyne  detector  when  shot 
noise  from  the  local  oscillator  is  the  dominant  noise  source  is  given  by  the 
relation 
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where  p  is  the  heterodyne  efficiency,  PR  is  the  signal  power  incident  upon 
the  detector,  hv  is  the  energy  of  a  signal  photon  with  frequency  v,  and  B  is 
the  larger  of  the  bandwidth  of  the  electrical  circuit  and  the  laser  bandwidth. 
Clearly,  the  transmitter  laser  should  possess  the  highest  possible  power  and 
the  narrowest  possible  bandwidth.  This  double  constraint  immediately  elimi¬ 
nates  many  of  the  possible  candidate  laser  systems.  Prominent  among  the 
remaining  candidates  are  visible  noble  gas  ion  lasers,  Nd3+:YAG  and  related 
solid  state  lasers,  HF  or  DF  lasers,  00  lasers,  and  CO^  lasers. 

The  visible  noble  gas  ion  lasers  are  very  inefficient,  requiring  kilowatts 
of  input  power  to  produce  a  few  watts  of  laser  output.  For  this  reason  they 
are  unsuitable  for  airborne  applications  where  the  available  electrical  power 
is  quite  limited.  Nd^+  lasers  can  be  quite  efficient  when  pumped  by  semi¬ 
conductor  diode  lasers.  However,  a  number  of  factors  related  to  their  short 
wavelength  mediate  against  their  use.  First,  since  hv  appears  in  the  denom¬ 
inator  of  Eq.  (1),  the  shorter  the  laser  wavelength  used  the  higher  the 
received  power  must  be  to  yield  the  same  SNR  (all  other  factors  being  equal) . 
Second,  efficient  heterodyne  detection  (maximum  n)  requires  that  the  signal 
and  local  oscillator  wavefronts  maintain  a  phase  relationship  which  is  con¬ 
stant  to  within  a  small  fraction  of  a  wavelength  over  the  entire  detector 
area.  Examination  of  Fig.  1  indicates  that  the  maximum  phase  mismatch 
between  the  wavefronts  of  two  collimated  beams  with  angular  separation  y 
incident  on  a  detector  of  linear  dimension  £  will  be 

A<t>  ■  y  2tt  *\#  2ir 
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Since  efficient  heterodyne  detection  requires  «  it,  the  shorter  the  incident 
wavelength  the  smaller  the  maximum  angular  separation  that  can  be  tolerated. 
Third,  the  absolute  frequency  stability  of  a  laser  becomes  harder  to  control 
as  the  wavelength  is  decreased.  This  is  complicated  by  the  large  gain  band- 
widths  of  Nd3+  lasers  which  further  increases  the  difficulty  of  stable  opera¬ 
tion.  Any  increase  in  the  bandwidth  B  results  in  a  decrease  in  SNR.  For 

3+ 

these  reasons  heterodyne  detection  of  Nd  and  other  visible  and  near  infra¬ 
red  lasers  is  a  lab  curiosity.  Longer  wavelength  lasers  are  better  when 
heterodyne  detection  is  employed  in  a  field  environment. 

Atmospheric  attenuation  is  another  factor  working  against  short  wave¬ 
length  lasers.  Fig.  2  shows  atmospheric  attenuation  data  obtained  for  four 
laser  wavelengths (C02  -  10.6  pm,  CO  -  4.73  pm,  DF  -  3.83  um,  and  Nd+3  :YAG  - 

4 

1.06  ym)  in  the  vicinity  of  Essen,  Germany.  In  winter,  the  attenuation  at 
10.6  ym  exceeds  2 dB/km  only  15%  of  the  time  while  the  attenuation  at  1.06  ym 
exceeds  2  dB/km  40%  of  the  time.  In  summer,  the  attenuation  at  10.6  ym 
exceeds  2  dB/km  only  8%  of  the  time  while  the  attenuation  at  1.06  ym  exceeds 
2  dB/km  18%  of  the  time.  Attenuation  at  wavelengths  below  1.06  ym  is 
characteristically  higher  still.  Although  attenuation  data  obtained  at 
other  sites  is  quantitatively  different,  qualitatively  the  attenuation  is 
always  much  less  at  10.6  v*n  faid  less  still  at  3.83  vm  or  4.73  urn)  than  at 
1.06  ym.  Because  a  small  increase  in  attenuation  drastically  reduces  the 
received  power  when  detection  ranges  of  several  kilometers  or  more  are 
involved,  and  in  view  of  the  arguments  presented  in  the  preceding  paragraph, 
we  conclude  that  Nd3+  lasers  are  unsuitable  for  use  in  the  IRAR  system. 
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Fig.  2.  Seasonal  atmospheric  attenuation  effects  in  Essen,  Germany  at  four 
laser  wavelengths:  00,  -  10.6  pm,  CO  -  4.73  pm,  DF  -  3.83  pm,  and  Nd3*:  YAG  - 
1.06  pm.  The  curves  plot  the  probability  that  the  attenuation  is  greater  than 
or  equal  to  the  value  given  on  the  abscissa. 
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Atmospheric  attenuation  arguments  favor  using  00  or  DF  lasers  over  002 
lasers.  Furthermore,  for  constant  transmitter  aperture  diameter,  the  angular 
resolution  of  the  radar  increases  with  decreasing  wavelength.  This  also 
favors  00  or  DF  lasers.  However,  00  lasers  operate  on  a  partial  population 
inversion  (an  inversion  between  specific  rotational  levels  of  different 
vibrational  states  without  an  inversion  of  the  total  population  of  the  two 
vibrational  states;  this  requires  a  rotational  temperature  which  is  signifi¬ 
cantly  lower  than  the  vibrational  temperature).  As  a  consequence,  efficient 
operation  of  CD  lasers  requires  cooling  to  liquid  nitrogen  temperatures . 

This  must  be  done  directly  in  a  conventional  electric  discharge  CD  laser  and 
would  be  prohibitive  in  an  airborne  application.  Cooling  is  achieved  by  gas- 
dynamic  processes  in  a  closed-cycle  GDL  00  laser.  However,  gasdynamic  lasers 
tend  to  be  quite  large  and  are  probably  unsuitable  for  use  on  a  tactical 
aircraft.  High-power  DF  lasers  are  typically  combustion  driven  systems  requir¬ 
ing  high  gas  flow  and  an  efficient  means  to  exhaust  the  DF  waste  product.  Such 
open-cycle  systems  are  impractical  for  use  on  tactical  aircraft.  Closed- 
cycle  electrically -pumped  DF  lasers  operate  very  similarly  to  CO  lasers  and 
thus  would  suffer  many  of  the  same  problems.  Therefore,  barring  significant 
technological  advances  in  CO  or  DF  laser  operation,  the  logical  choice  of 
laser  system  for  use  in  the  IRAR  MTI  subsystem  is  a  C02  laser. 

Longitudinal  electric  discharge  C02  lasers  have  long  been  proven  capable 
of  achieving  high  powers  (in  either  cw  or  repetitively -pulsed  modes) ,  narrow 
linewidths,  and  efficient  heterodyne  detection.  Furthermore ,  they  can  be 
made  sufficiently  compact  and  efficient  to  be  readily  adaptable  to  airborne 
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use. 3  In  an  MTI  application  the  question  immediately  arises  as  to  whether  cw 
or  repetitively-pulsed  operation  is  optimum.  At  C02  laser  frequencies 
(>  3  x  1013  Hz)  the  Doppler  shift  corresponding  to  2  kmph  target  motion  is 
only  106  kHz.  To  be  able  to  achieve  this  frequency  resolution  the  effective 
laser  bandwidth  must  be  less  than  106  kHz  and  consequently  the  laser  radiation 
must  have  a  duration  which  is  in  excess  of  12  usee.  For  a  cw  laser  the 
effective  duration  is  given  by  the  single  pixel  dwell  time.  As  will  be  shown 
later  this  can  be  in  excess  of  100  ysec.  For  dwell  times  in  excess  of  the 
minimum  allowable  duration,  repetitively -pulsed  operation  can,  in  principle, 
result  in  increased  peak  powers  with  consequent  increased  signal -to -noise 
ratios  without  having  the  pulse  bandwidth  exceed  the  desired  Doppler  resolu¬ 
tion.  However,  in  practice,  it  is  exceedingly  difficult  to  construct  compact, 
efficient,  high  duty  cycle,  repetitively -pulsed  C02  lasers  with  the  desired 
bandwidth  and  pulse  repetition  frequency  (PRF)  characteristics.  Furthermore, 
as  will  also  be  shown  later,  much  of  the  increase  in  SNR  resulting  from  pulsed 
operation  can  also  be  realized  by  suitable  post-detection  integration  of  the 
cw  laser  signal.  Thus,  a  cw  C02  laser  is  the  optimum  choice  for  use  in  the 
IRAR  MTI  subsystem. 

The  choice  of  appropriate  detector  material  is  fairly  straightforward. 
When  high  frequency  heterodyne  detection  is  involved,  photovoltaic  semicon¬ 
ductor  detectors  are  far  superior  to  photoconduct ive  detectors.  Furthermore, 
detector  cooling  requirements  should  be  kept  to  a  minimum.  At  10.6  nn  the 
only  materials  available  as  photovoltaic  detectors  which  can  be  operated  at 
temperatures  as  high  as  77°K  (liquid  nitrogen)  are  Hg^^Cd^Te  and 
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pk(l-x)SnxTe*  These  two  materials  have  very  similar  properties,  both  can  be 
fabricated  into  arrays,  and  both  make  excellent  heterodyne  detectors  at 
10.6  pm  with  high  quantum  efficiencies.  Either  would  be  suitable  for  use  in 
the  IRAR  MTI  subsystem. 

3.  SIZE  AND  SHAPE  OF  THE  HETERODYNE  DETECTOR  ARRAY 

Individual  HgCdTe  detectors  can  be  made  quite  small  and  can  be  readily 
incorporated  into  large  one-  and  two-dimensional  detector  arrays.  A  number 
of  factors  influence  whether  a  single  detector  will  suffice  for  use  in  the 
IRAR  MTI  subsystem  or  whether  a  detector  array  is  needed.  If  an  array  is 
required,  still  other  factors  will  determine  the  optimum  number  of  detectors 
and  the  geometry  of  their  arrangement  in  an  array.  In  this  section  we  will 
examine  the  requirements  iiqposed  on  the  detector  system  by  the  IRAR  MTI 
mission. 

Consider  the  system  shown  in  Fig.  3.  The  linear  size  of  a  laser  beam  of 
wavelength  X  transmitted  by  an  aperture  of  diameter  Dp  onto  a  target  at  range 
R  (assuming  diffraction-limited  divergence  only)  is  given  by 

L-p  %  Rap  R;  RX/Dtp  (3) 

otp  denotes  the  angular  resolution  of  the  transmitter  aperture .  Approximately 
equal  signs  are  used  in  Eq.  (3)  and  succeeding  equations  because  the  exact 
numerical  coefficients  (all  of  order  unity)  depend  on  the  exact  intensity 
profiles  of  the  transmitted  beams.  The  phase  correlation  distance,  L^, 
observed  at  a  distant  receiver,  of  the  radiation  scattered  from  a  diffusely 
reflecting  surface  is  related  to  the  size  of  the  illuminated  region  by6 
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Lc  %  XR/Lj.  *  Up 


(4) 


If  the  receiver  aperture  (of  diameter  Dr)  is  larger  than  Lc,  the  lack  of  phase 
correlation  across  the  receiver  aperture  prohibits  efficient  heterodyne  detec¬ 
tion  with  a  single  detector  of  all  of  the  radiation  scattered  from  the  spot  of 
size  Ly. 

Consider  next  a  small  portion  of  the  illuminated  spot.  The  phase  correla¬ 
tion  distance  at  the  receiver  corresponding  to  a  portion  of  diameter  Li,  is 


LI  %  AR/LJ.  (5) 

Thus  a  single  detector  matched  to  a  receiver  of  diameter  Dr  ^  can  effi¬ 
ciently  heterodyne  the  radiation  from  a  portion  of  the  scene  of  size  Lj,. 

The  angular  separation  of  adjacent  portions  in  the  scene  of  size  LJ.  is  given 
by 

Or  £  LJ.  /R  £  A/Dr  (6) 

or  just  the  angular  resolution  of  an  aperture  of  diameter  Dr.  Thus,  if 
we  consider  a  receiver  of  diameter  Dr  and  place  an  array  of  closely -spaced 
detectors  at  the  focal  plane  of  a  lens  imaging  the  receiver  aperture  such 
that  each  detector  looks  at  a  different  portion  of  the  illuninated  scene  of 
size  LJ. ,  an  efficient  heterodyne  output  can  be  obtained  from  each.  The 
number  of  detectors  required  to  completely  cover  the  illuminated  scene  is 
given  roughly  by  the  ratio  of  the  illuminated  area  to  the  area  viewed  by  a 


single  detector,  i.e.. 


If  it  is  desired  to  obtain  a  single  value  from  the  whole  illuminated  spot,  the 
outputs  of  the  N  detectors  must  be  incoherently  added.  Coherent  addition 
of  the  outputs  is  not  viable  because  the  heterodyne  signals  retain  information 
on  the  optical  phase.  The  lack  of  phase  correlation  across  the  aperture  that 
prevents  efficient  heterodyne  detection  with  a  single  detector  results  in  a 
lack  of  phase  correlation  between  the  N  detectors  which  prevents  coherent 
addition  of  the  outputs. 

The  preceding  results  were  derived  for  target  plane  spot  sizes  produced 
solely  by  diffraction  from  the  apertures.  However,  they  are  equally  valid 
for  spot  sizes  produced  by  artificially  diverging  (with  a  negative  lens)  the 
beam  coming  from  a  larger  aperture.  As  the  aperture  diameters  are  no  longer 
the  factors  determining  the  divergences,  the  forms  of  the  preceding  expres¬ 
and  beam  divergences  (otp, 

alone  should  be  employed  for  calculational  purposes . 

Now  let  us  consider  an  aircraft  in  level  flight  with  velocity  V  and  alti¬ 
tude  H  projecting  a  laser  beam  onto  the  ground  at  slant  range  R  (see  Fig.  4). 
Let  the  transmitter  have  an  aperture  of  diameter  Up  and  the  transmitted  beam 
have  (lens -induced)  divergences  and  cy  *  Ny,  where  ctp  ■  X/Dtp, 

Nx  and  Ny  are  constants  (>  1) ,  x  is  the  direction  perpendicular  to  the  plane 
defined  by  the  altitude  and  slant  range  vectors,  and  y  is  the  direction 


sions  which  involve  the  spot  sizes  (Lp,  Lp 
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Fig.  4.  Geometry  encountered  in  the  tactical  aircraft  application  of  the 
IRAR  system. 


perpendicular  to  both  the  x-direction  and  the  slant  range  vector.  The  physi¬ 
cal  size  of  the  beam  when  projected  onto  the  ground  is 

Lx  =  2axR  =  2NxarR  (g) 

for  the  x-direction  and 

Ly  “  2°>yR/sinft  *  ZNyOjRfR/H)  (9) 

for  the  y- direction. 

The  time  required  for  the  aircraft  to  move  forward  a  distance  equal  to 

Ly  is 

t  -  Ly/v  =  2NyapR2/HV  (10) 

If  we  assune  that  the  beam  is  swept  along  the  ground  in  the  x-direction  by  a 
scanner,  the  time  required  to  complete  a  scan  with  a  field-of-view  ex  is 

tx  -  Sj/Bw  (ID 

where  Jj  is  the  angular  scan  rate  and  8  is  the  retrace  efficiency  (scan  time 
divided  by  scan  time  plus  retrace  time) .  In  a  search  mode  it  is  desirable 
to  illuninate  every  portion  of  the  ground  (within  the  x-direction  field-of- 
view)  at  least  once.  This  is  accomplished  when  t  >  t  ,  that  is,  when 

&  >  e^/ZSNyOpR2  (12) 

Typical  scanners  are  limited  to  modest  scan  rates,  thus  for  some  values  of 
the  parameters  9X,  8,  H,  V,  Op,  R,  and  u  it  may  be  necessary  to  have  Ky  »  1 
in  order  to  satisfy  Eq.  (12). 
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Assume  that  the  illuminated  spot  is  viewed  by  a  receiver  with  aperture 
D^.  A  detector  array  consisting  of  nx  x  detectors,  where  nx  and  iiy  are 
the  minimum- valued  integers  such  that 

%  1  lWnr  <13> 

and 

ny  1  NyVDT>  (14) 

is  required  to  efficiently  heterodyne  detect  the  radiation  scattered  by  the 
spot.  If  all  other  parameters  are  held  constant  the  signal -to -noise  ratio 
of  the  output  of  each  detector  is  inversely  proportional  to  the  number  of 
detectors 

SNR  =  8’Ayy  (15) 

because  each  detector  looks  at  only  a  fraction  of  the  illuminated  spot  and 
thereby  receives  only  the  laser  radiation  scattered  from  that  fraction.  6' 
contains  the  dependences  on  laser  power,  optical  efficiency,  collection  solid 
angle,  atmospheric  attenuation,  and  target  reflectivity. 

Targets  of  interest  typically  encompass  a  number  of  receiver  resolution 
elements  or  pixels.  To  allow  maximum  operational  utility  of  the  system  in 
bad  weather,  the  system  must  have  the  highest  possible  SNR.  Incoherent 
integration  (addition)  of  the  outputs  from  a  number  of  pixels  (all  of  which 
lie  on  the  target)  yields  an  improvement  in  the  SNR.  Specifically,  incoher¬ 
ent  integration  of  the  outputs  of  n'x  detectors  along  the  scan  direction  in 
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each  row  of  the  array  yields  an  improvement  Gx 

1  <  Cx  <  n*x  (nx  >  n'x  >  1)  (16) 

while  incoherent  integration  of  the  outputs  of  n'y  detectors  in  each  column  of 
the  array  yields  an  inprovement  Gy 

1  £  Gy  <  n'y  (ny  >  n'y  >  1)  (17) 

Furthermore ,  each  detector  looks  at  any  given  pixel  for  one  dwell  time 

tD  *  X/Dj^d)  (18) 

but  is  on  an  object  n'x  pixels  wide  for  n'x  dwell  times.  Consequently,  tem¬ 
poral  integration  of  the  signal  from  each  detector  over  n'x  dwell  times  gives 
a  SNR  improvement  Gt 

1  <  Gt  <  n'x  (19) 

Thus  the  final  SNR  which  can  be  obtained  for  a  two-dimensional  array  with  an 

effective  resolution  of  n*  x  n'  pixels  is 

x  y 

SNR2D  *  S'CxGyVVy  (20) 

In  the  special  case  of  a  one -dimensional  array  (r^  »  1),  no  spatial 
integration  can  be  performed  in  the  x-direction  and  Gx  *■  1.  As  with  the 
two-dimensional  array  spatial  integration  in  the  y-direction  yields  a  gain 
given  by  Eq.  (17)  while  tenporal  integration  yields  a  gain  given  by  Eq.  (19). 
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Thus,  assuming  that  the  transmitted  beam  is  shaped  to  match  a  one -dimensional 
array  according  to  Eqs.  (13)  and  (14),  the  final  SNR  vhich  can  be  obtained 
for  a  one -dimensional  array  with  an  effective  resolution  of  n*  x  n*  pixels 

x  y 

is 

SNR1D  =  S'GyGt/ny  =  (njGx)  SNR2D  (21) 

Since  from  Eq.  (16)  Gx  <  n'  <_  nx,  we  find  that  SNRjD  >  SNR2D.  Each  detector 
in  the  array  must  have  its  own  Doppler  analyzer.  Therefore,  the  two-dimen¬ 
sional  array  requires  nx  times  as  many  analyzers  as  the  one -dimensional  array. 
As  a  result  we  find  that  a  one -dimensional  array  is  superior  to  a  two-dimen¬ 
sional  array  in  a  line  scan  MTI  application  for  reasons  of  both  system 
simplicity  and  increased  signal -to-noise  ratio. 

Efficient  heterodyne  detection  with  a  one -dimensional  array  requires 
that  the  spot  size  viewed  by  the  receiver  aperture  be  larger  than  the  trans¬ 
verse  size  of  the  illuminated  spot.  Thus  we  find  that  we  must  have 

W“r  i  1  (22) 

Since  the  signal  received  by  a  detector  increases  as  the  square  of  the 
receiver  diameter,  it  is  clear  that  we  want  to  make  as  large  as  possible. 
This  coupled  with  the  fact  that  N  >  1  implies  that  the  optimum  situation 

A  *** 

arises  when 

*  Dtp  and  Nx  ■  1 

i.e.,  when  there  is  no  lens -induced  divergence  in  the  x-direction  and  the 
receiver  and  transmitter  apertures  are  of  equal  size.  Given  this  situation 
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both  system  simplicity  and  system  conpactness  then  dictate  that  the  trans¬ 
mitter  and  receiver  should  share  the  same  aperture. 

The  same  arguments  that  favor  a  one -dimensional  array  over  a  two-dimen¬ 
sional  array,  also  favor  a  single  detector  over  a  one -dimensional  array. 
However,  as  shown  earlier,  complete  coverage  of  the  scene  may  require  Ny »  1. 
If  Ny  »  1  and  =  Up,  then  a  single  detector  is  insufficient  and  a  one¬ 
dimensional  array  must  be  used. 

The  complete  expression  governing  the  signal -to -noise  ratio  obtained  for 


any  individual  detector  in  a  one -dimensional  array  is 


.  1 


SNR 


?t  ^  Pone"2kR 

— r  e p ne 

ATk*- 


NjjhvB  4R 


(23) 


where  PT  is  the  transmitted  laser  power,  Np  is  the  total  number  of  detectors 

in  the  array,  e  is  the  efficiency  of  the  optical  system,  p  is  the  target 
-2kR 

reflectivity,  e  is  the  round  trip  atmospheric  attenuation,  and  the  remain¬ 
ing  parameters  have  their  previously-defined  meanings.  Because  we  require 
complete  coverage  of  the  ground  within  the  horizontal  field-of-view  in  the 
line  scan  mode  of  operation,  the  receiver  aperture,  the  number  of  detectors, 
and  the  required  bandwidth  are  all  interdependent.  That  is,  if  the  minimum 
number  of  detectors  required  to  cover  the  conplete  scene  for  a  given  and 
u  are  used,  changing  any  one  of  Np,  D^,  &,  or  B  will  cause  changes  in  the 
others.  If  we  separate  out  all  of  the  interdependent  factors,  we  may  write 
Eq.  (23)  as 

2 

SNR  -  (24) 
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If  the  system  bandwidth  is  limited  by  the  dwell  time  of  the  laser  beam  on  a 
single  pixel,  then  from  Eqs.  (12)  and  (18)  we  have 

B  *  Bd  =  l/tD  =  =  tI)D^/2X 

*  0xHTOr/26XNdoR2  =  0XHVD2R/28A2R2ND  (25) 

Substituting  Eq.  (25)  into  Eq.  (24)  we  find 

2  2 

SNR  -  B'  =  28'BA  R  /e  HV 

%% 

_  PT  BA2  -2kR  _  CMD 

“  Hv  epne  =  SNRo  (26) 

Thus,  the  signal -to-noise  ratio  obtained  when  the  MTI  system  is  dwell  time 
band -width- limited  is  independent  of  the  receiver  aperture  or  the  scan  rate. 

It  is  determined  by  the  laser  power,  the  mission  requirements  (H,  V,  0x,  R, 
and  k) ,  and  the  system  efficiency.  This  conclusion  does  not  contradict  any 
of  our  earlier  statements,  because  earlier  we  allowed  the  possibility  that 
f  Up,  a  possibility  that  was  subsequently  shown  to  be  non-optimum. 

If  the  bandwidth  of  the  system  is  not  dwell  time-limited  but  instead 
is  limited  by  laser  stability  (or  minimum  obtainable  Doppler  filter  bandwidth) 
then 

B  =  Bl 
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In  this  case  Eq.  (24)  becomes 


s*  '  fr  ’  ®d  i  \  «7> 

nd\ 

That  is,  the  signal -to-noise  ratio  in  the  laser  bandwidth- limited  regime  is 
lower  than  that  which  would  be  obtained  for  a  dwell -time  bandwidth- limited 
system  with  identical  and  <J.  The  behavior  of  Eqs.  (26)  and  (27)  as  a 
function  of  Bp/Bj^  is  summarized  in  Fig.  5.  For  Bp  <  B.,  the  SNR  is  laser 
bandwidth- limited  and  Eq.  (27)  is  valid,  while  for  B^  >_  B^^  the  SNR  is  dwell 
bandwidth-limited  and  obeys  Eq.  (26).  It  should  be  noted  that  this  behavior 
is  characteristic  of  the  optimized  MTI  line  scan  mode  of  operation  and  would 
not  necessarily  be  valid  for  any  other  mode  of  IRAR  operation. 

As  mentioned  earlier,  fluctuations  in  the  effective  target  reflectivity 
(due  to  speckle  or  atmospheric  turbulence)  will  reduce  the  MTI  detection 
probability.  Thus  it  is  desirable  to  integrate  the  returns  from  a  nunber  of 
pixels  (all  of  which  lie  on  the  target)  to  increase  the  detection  probability. 
This,  of  course,  requires  that 

2apR  -  2XR/Qj.  «  L  (28) 

where  L  is  a  characteristic  target  dimension.  For  targets  of  2  m  character¬ 
istic  dimension  located  at  ranges  of  the  order  of  3  km,  Eq.  (28)  is  satisfied 
for  aperture  diameters  of  the  order  of  10  cm  or  larger. 

Other  factors  also  affect  the  choice  of  aperture  size.  First,  the  fact 
that  the  entire  system  must  be  capable  of  being  carried  by  a  tactical 
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Fig.  5.  Summary  of  bandwidth  effects  on  the  signal -to -noise  ratio  of  an 
MTI  system.  Bo  is  the  dwell  bandwidth,  Bl  is  the  laser  bandwidth,  and  SNRq 
is  the  dwell  bandwidth-limited  signal-to-noise  ratio. 


aircraft  will  in  practice  limit  the  maxi  man  aperture  that  can  be  used.  In 
addition,  the  larger  the  aperture,  the  larger  the  scanner  mirror  that  is 
required  and  the  slower  the  maximun  scan  rate  that  can  be  achieved  in  practice. 
Slower  scan  rates  coqpled  with  smaller  resolution  angles  resulting  from 
increased  aperture  sizes  nean  that  more  detectors  are  required  for  adequate 
coverage  of  the  scene.  This  in  turn  means  an  increase  in  the  size  and  com¬ 
plexity  of  the  NTTI  electronic  processing  system.  Consideration  of  all  of 
the  above  factors  indicate  that  the  aperture  of  the  IRAR  MTI  system  should 
fall  somewhere  in  the  range  of  10-30  cm  diameter.  Finally,  it  is  highly 
desirable  that  the  IRAR  MTI  subsystem  be  compatiable  with  the  IRAR  imager 
subsystem.  Since  the  imager  aperture  has  been  sized  at  approximately  10  cm 
for  angular  resolution  reasons,  we  will  henceforth  assume  this  value  for  the 
MTI  system  as  well. 

Given  all  of  the  preceding  information  we  are  now  in  a  position  to  esti¬ 
mate  the  size  of  detector  array  that  is  required  by  the  IRAR  MTI  system. 
Maximum  scan  rates  for  a  10  cm  diameter  full  aperture  scanner  are  probably 
limited  to  values  of  the  order  of  1  rad/sec  with  70%  retrace  efficiency.  If 
we  consider  an  aircraft  flying  at  150  m/sec  Mach  0.5)  at  50  m  altitude 
(values  conpatible  with  a  close  air  support  mission)  and  searching  a  1  radian 
field-of-view  at  3  km  slant  range,  we  find  from  Eq.  (12) 

It  is  interesting  to  note  that  a  13-element  one -dimensional  array  is  very 
close  to  the  12 -element  array  proposed  for  the  IRAR  imager  system.  For  the 
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above  system  parameters,  from  Eq.  (18)  we  find  a  dwell  time  of 


tD  =  100  usee. 


with  a  corresponding  bandwidth  of  10  kHz. 

4.  SIGNAL-TO-NOISE  RATIO  AND  CLUTTER-TO-NOISE  RATIO  CONSIDERATIONS 

The  ability  of  any  MTI  system  to  detect  a  moving  target  is  profoundly 
influenced  by  system  noise.  Because  the  optimum  technique  for  detecting 
the  presence  of  a  signal  in  noise  is  a  simple  threshold  test,  random  noise 
spikes  can  exceed  the  threshold  level  and  give  false  alarms.  If  the  thresh¬ 
old  is  raised  high  enough  to  limit  false  alarms  to  an  acceptable  rate,  low- 
level  signals  may  not  be  strong  enough  to  exceed  the  threshold  and  the 
detection  probability  will  be  reduced.  The  interrelationships  of  the  detec¬ 
tion  probability  Pp,  the  false  alarm  probability  Pp,  and  the  signal -to-noise 
ratio  SNR  are  summarized  in  receiver  operating  characteristic  (ROC)  curves. 

For  targets  obeying  speckle  statistics  (i.e.,  diffusely  reflecting  tar¬ 
gets),  Pp,  Pp,  and  SNR  are  related  by  the  expression® 


p  .  p  (1  ♦  SNR) 
D  "  *F 


(29) 


ROC  curves  appropriate  to  Eq.  (29)  are  plotted  in  Fig.  6  for  several  values 
of  Pp. 

For  an  MTI  system  with  Np  detectors  and  N^  Doppler  channels  of  bandwidth 
B^.  per  detector,  the  false  alarm  probability  required  in  each  channel  is 
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g.  6.  Receiver  operating  characteristic  curves  for  speckle  targets. 
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1/N^c  times  the  false  alarm  probability  desired  for  the  system  as  a  whole. 
Thus, 

pf  ’  “TO  (io) 

where  Tp  is  the  required  mean  time  between  false  alarms  for  the  system  as  a 
whole.  If  the  Doppler  signal  returned  from  a  single  pixel  has  a  bandwidth 
which  is  less  than  the  Doppler  channel  bandwidth,  then  the  SNR  calculated 
from  Eq.  (23)  may  be  conpared  directly  with  the  abscissa  of  Fig.  6.  This  con¬ 
dition  should  hold  for  a  properly -designed  <X>2  laser  MTI  system.  If  it  is 
not  valid,  the  result  of  Eq.  (23)  must  be  reduced  by  the  ratio  of  the  band- 
widths.  This  occurs  because  any  one  Doppler  channel  will  be  capable  of  pass¬ 
ing  only  a  fraction  of  the  total  signal. 

Typical  parameters  for  the  IRAR  system  are  sunmarized  in  Table  I.  Using 
these  values  and  ass  lining  an  average  target  reflectivity  p  of  0.1  and  an 
atmospheric  attenuation  of  2  dB/km  (this  or  a  lesser  value  will  occur  65%  of 
the  time  in  Central  Europe1,4),  Eq.  (23)  yields 

SNR  15.5  -  11.9  dB 

For  false  alarm  times  Tp  of  the  order  of  15  minutes,  we  obtain  from  Eq.  (30) 

Pc  %  10"11 
F 

Consulting  Fig.  6  we  find  that  the  system  signal-to-noise  ratio  is  woefully 
inadequate  to  achieve  detection  probabilities  greater  than  901  with  Pp  =  10  11 
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TABLE  I 


TYPICAL  PARAMETERS  FOR  THE  IRAR  SYSTEM 


Laser  power  -  P 
Photon  energy  -  hv 
Number  of  detectors  -  Np 
Number  of  Doppler  channels 
Bandwidth  per  channel  - 
Dwell  bandwidth  -  Bp 
Aperture  diameter  -  D^ 
Maximum  range  -  R 
Optical  efficiency  -  e 
Heterodyne  Efficiency  -  n 


10  W  cw 
1.87  x  10"20 
12 

-  Nc  100 

100  kHz 
10  kHz 
10  cm 
3  km 
0.1 


0.5 


For  speckle  targets  the  receiver  operating  characteristics  at  high  Pp 
values  are  dominated  by  the  probability  that  speckle  will  cause  a  deep  fade 
in  the  received  signal  with  a  resultant  failure  to  exceed  the  threshold  level. 
If  several  independent  pixels,  all  of  which  lie  on  the  target,  are  examined, 
the  probability  that  at  least  one  of  the  returns  will  exceed  the  threshold  is 
greatly  increased.  If  N  returns  are  examined,  the  detection  probability 
becomes 

Pd(N)  =  1  -  C1-Pf(1  +  SNR)  V  (31) 

The  dashed  curve  in  Fig.  7  shows  the  number  of  pixels  which  must  be  examined 

-12 

to  yield  P^  =  0.99  and  Pp  =10  as  a  function  of  the  single-pixel  signal- 

to -noise  ratio.  For  SNR  =  11.9  dB,  at  least  22  pixels  must  lie  on  the  target. 

Since  a  typical  tactical  target  at  3  km  range  will  occupy  50-100  pixels,  some 

degree  of  multiple -pixel  processing  is  certainly  allowable. 

However,  simply  setting  the  threshold  to  give  the  desired  single-pixel 

Pp  value  and  hoping  that  the  number  of  pixels  on  the  target  is  large  enough 

so  that  at  least  one  return  exceeds  the  threshold  is  by  no  means  optimal.  A 

much  more  efficient  method  is  to  incoherently  add  the  returns  from  N  pixels 

and  compare  the  integrated  signal  with  a  threshold  which  is  a  function  of  the 

nunber  of  pixels  integrated.  The  Pp/Pp  effects  of  pixel  integration  can  be 

g 

mathematically  evaluated  using  moment  generating  function  techniques.  The 
solid  curve  in  Fig.  7  shows  the  number  of  pixels  which  must  be  integrated  to 
achieve  PD  =  0.99  and  Pp  =  lO’1^  as  a  function  of  the  single-pixel  SNR.10 
For  SNR  *  12  dB,  pixel  integration  requires  only  7  pixels  as  opposed  to  22 
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Fig.  7.  Speckle  target  multiple-pixel  processing  inprovement  curves  for 
Pp  ■  0.99  and  Pp  ■  10‘i2.  The  dashed  curve  gives  the  number  of  target  returns 
which  must  be  examined  as  a  function  of  the  single-pixel  SNR  if  the  threshold 
is  set  to  yield  the  desired  Pp  on  a  single-pixel  basis  and  no  integration  is 
performed.  The  solid  curve  gives  the  nunber  of  pixels  which  must  be  integrated 
if  the  threshold  is  set  to  be  optimum  for  a  signal  consisting  of  N  incoherently 
added  target  returns. 


pixels  in  the  single-pixel  threshold  approach.  It  should  be  noted  that  the 
curves  in  Fig.  7  are  not  strong  functions  of  either  Pp  or  Pp. 

The  preceding  results  are  valid  for  speckle  (diffusely-reflecting)  tar¬ 
gets.  Some  objects,  however,  exhibit  specular  reflections  called  glints. 
Glints  do  not  fluctuate  like  speckle  targets  and  therefore  will  have  different 

g 

Pp/Pp  behavior.  The  detection  probability  for  a  glint  target  is  given  by 

PD  =  Qp[(2  SNR) 1/2 ,  (-2  £nPp)1/2]  (32) 

Q 

where  Qp[u,  v]  is  Marcum's  Q  function.  ROC  curves  appropriate  to  Eq.  (32) 
are  plotted  in  Fig.  8 Comparison  of  Figs.  6  and  8  indicates  that  much 
lower  signal -to -noise  ratios  are  required  to  detect  glint  targets  than  to 
detect  speckle  targets  at  high  Pp  values.  This  suggests  that  the  presence 
of  target  glints  may  reduce  the  need  for  multiple -pixel  processing.  However, 
as  we  shall  see,  this  is  only  partially  true. 

If  the  threshold  is  chosen  to  produce  the  desired  single-pixel  Pp  and 
N  pixels  (all  containing  a  glint  target)  are  examined  without  integration, 
the  effective  detection  probability  is  given  by 

Pp(N)  =  1  -  (1-Q1[(2  SNR)1/2,  (-2  £,nPp)1/2])N  (33) 

Eq.  (33)  has  been  used  to  evaluate  the  number  of  pixels  required  to  obtain 
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Pp  =  0.99  with  Pp  ■  10  for  various  SNR  values.  These  results  are  summar¬ 
ized  in  the  dashed  curve  of  Fig.  9.  As  with  speckle  targets  pixel  integration 
produces  much  better  results  for  glint  targets  than  does  the  single -pixel 
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INITIAL  SINGLE-PIXEL  SNR  (dB) 

Fig.  9.  Glint  target  multiple-pixel  processing  improvement  curves  for  P™  * 
0.99  and  Pp‘12.  The  dashed  curve  represents  results  obtained  for  the  single- 
pixel  threshold  approach.  The  solid  curve  represents  results  obtained  for  the 
j  pixel  integration  approach. 

I 
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threshold  approach.  Using  moment  generating  function  techniques  the  solid 
curve  of  Fig.  9  was  obtained. ^  Close  comparison  of  Figs.  7  and  9  shows  that 
for  glint  targets  multiple -pixel  processing  is  not  required  for  SNR  >  16  dB 
while  for  speckle  targets  this  condition  is  only  satisfied  for  SNR  >  33  dB. 
However,  for  SNR  <  10  dB  conparable  nunbers  of  pixels  must  be  processed  for 
either  glint  or  speckle  targets.  Thus,  if  the  target  glints  produce  returns 
which  are  much  stronger  than  the  expected  speckle  returns  (>  10-15  dB) 
pixel  integration  will  not  be  required,  while  if  the  glint  returns  are  only 
of  the  same  average  strength  as  the  expected  speckle  returns,  pixel  integra¬ 
tion  will  be  a  necessity.  Since  the  presence  of  strong  target  glints  cannot 
be  counted  on,  some  degree  of  pixel  integration  must  be  incorporated  into  the 
MTI  system. 

If  the  Doppler  channel  bandwidth  Bc  is  significantly  larger  than  the 
dwell  bandwidth  Bp,  an  additional  improvement  in  signal-to-noise  ratio  may  be 
realized  by  performing  a  postdetection  integration  operation  on  each  channel 
before  pixel  integration.  This  is  possible  because,  in  heterodyne  detection, 
target  speckle  will  cause  the  signal  to  vary  on  time  scales  of  the  order  of 
1/Bp  while  the  noise  contribution  will  vary  on  time  scales  of  the  order  of 
1/Bj,.  Thus,  in  general,  nt  independent  signal -plus -noise  measurements  may  be 
made  during  the  time  it  takes  for  the  laser  beam  to  scan  across  a  single  pixel, 
where 


-  Bc/Bp 


(34) 
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The  inprovement  in  signal -to-noise  ratio  which  may  be  obtained  by  integrating 

nt  measurements  after  envelope  detection  of  the  IF  signal  is  plotted  in 

Fig.  10.  This  curve  is  derived  from  data  in  Ref.  11  for  Pp  *  0.99  and  Pp  * 
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10  .  Again  these  results  are  quite  weak  functions  of  Pp  and  Pp. 

Depending  on  the  magnitude  of  nt  the  additional  SNR  gained  by  postdetec¬ 
tion  integration  may  significantly  reduce  the  nunber  of  pixels  which  must  be 
averaged.  In  our  nunerical  example  of  Table  I,  Bq/Bq  *  10,  which  yields  a 
SNR  inprovement  of  approximately  9  dB.  This  would  reduce  the  speckle  averag¬ 
ing  requirement  from  7  to  2  pixels.  Since  postdetection  integration  still 
does  not  eliminate  the  need  for  pixel  integration  in  most  cases,  the  optimum 
approach  is  to  use  postdetection  integration  and  integration  of  as  many  pixels 
as  possible  while  insuring  that  the  number  of  pixels  integrated  is  less  than 
the  number  of  pixels  on  the  target.  For  tactical  targets  at  ranges  of  inter¬ 
est  and  allowing  for  different  target  orientations,  the  maximum  nunber  of 
pixels  which  may  be  integrated  is  approximately  25-30.  If  in  the  numerical 
example  30  pixels  were  integrated  instead  of  2,  an  additional  15  dB  of  SNR 
margin  is  obtained.  This  margin  would  allow  the  system  to  operate  out  to 
3  km  even  if  the  atmospheric  attenuation  were  as  high  as  4.5  dB/km.  Thus 
maximizing  the  nunber  of  pixels  integrated  maximizes  the  operational  utility 
of  the  MTI  system. 

It  should  be  noted  that  postdetection  integration  of  n^  =  n^  nt 
measurements  corresponds  to  pixel  integration  of  n^  pixels  for  each  detector 
of  a  one -dimensional  array,  in  addition  to  postdetection  integration  of  n^ 
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Fig.  10.  Signal-to-noise  ratio  improvement  obtained  by  postdetection  integra 
tion  of  n*  independent  signal -plus -noise  measurements  for  Pn  -  0.99  and 

PF  .  10-1 1. 
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measurements •  It  also  corresponds  to  the  temporal  integration  gain  Gt 
mentioned  in  the  preceding  section. 

Ground  clutter  can  have  a  significant  effect  on  the  performance  of  any 
Doppler  MTI  system.  Windblown  trees,  shrubbery,  and  water  waves  can  produce 
radar  returns  with  non-zero  Doppler  shifts  which  might  be  mistaken  for  slowly- 
moving  targets.  Furthermore,  because  any  filter  technique  has  only  a  finite 
rejection  for  signals  outside  its  bandpass,  filter  leakthrough  of  zero  velocity 
signals  can  affect  the  effective  noise  level  and  thus  the  effective  false 
alarm  probability. 

The  analysis  of  clutter  effects  is  generally  very  complicated.  However, 
by  making  certain  assumptions  the  problem  can  be  simplified  considerably. 

First,  we  assume  that  clutter  and  target  returns  are  mutually  exclusive.  That 
is,  the  IF  signal  received  from  each  pixel  consists  of  either  clutter -plus - 
noise  or  target -plus -noise,  but  not  clutter-plus -target -plus -noise.  In  a 
high  resolution  radar  (capable  of  imaging  objects),  this  is  a  good  approxima¬ 
tion.  Second,  we  assume  that  both  the  target  and  the  clutter  returns  obey 
speckle  statistics.  This  represents  a  worst  (although  highly  probable)  case 
from  the  target  standpoint,  and  a  typical  case  from  the  clutter  standpoint. 

Of  all  natural  clutter  sources  only  water  waves  are  expected  to  produce  moving 
glint  returns. 

For  a  speckle  target  the  detection  probability  P^p  is  related  to  the 
false  alarm  probability  (caused  by  noise  alone)  Pp  and  the  target  signal -to- 
noise  ratio  SNR  by  Eq.  (29) .  For  radar  returns  from  clutter  obeying  speckle 
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statistics  and  having  clutter  signal -to-noise  ratio  CNR,  the  detection  proba¬ 


bility  Ppj,  is  given  by 


p  .  p  (1  +  CNR) 
*DC  *F 


However,  any  detected  clutter  return  is  essentially  a  false  alarm.  Therefore, 
in  a  realistic  (i.e.,  cluttered)  scenario  we  may  define  an  effective  false 


alarm  probability 


P*F  =  PDC 


Combining  Eqs.  (29),  (35),  and  (36)  we  find  for  SNR  »  1 


f  1  .  CNR-, 

p,  _  p  (1  +  SNR)/ (1  +  QJR)  p  lSRR:+W 
*F  “  MJT  DfT 


The  behavior  of  P£  versus  the  clutter-to-signal  ratio  CNR/ SNR  is  plotted  in 

Fig.  11  for  target  detection  probability  P^  *  0.99  and  several  values  of  the 

noise-only  false  alaim  probability  Pp.  Fig.  11  indicates  that  clutter  returns 

must  be  kept  at  least  SO  dB  below  the  target  returns  or  the  effective  false 

alarm  probability  will  be  significantly  increased  above  its  noise-only  value. 

The  preceding  results  are  valid  for  single-pixel  detection.  However, 

we  have  already  seen  that  pixel  integration  is  required  to  solve  the  target 

detection  problem.  Using  moment  generating  function  techniques  the  effect 

of  pixel  integration  on  the  Pp  versus  CNR/ SNR  curves  can  also  be  evaluated.10 

In  Fig.  12  we  conpare  curves  for  integration  of  1,  10,  and  50  pixels  with 
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PDT  *  0  .99  and  Pp  ■  10  .  Even  with  pixel  integration  clutter  must  still 
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Fig.  11.  Effects  of  varying  the  clutter-to-signal  level  on  the  effective 
false  alarm  probability.  Curves  a^e  plotted  for  various  noise-only  false  alarm 
probabilities  (Pc)  and  assume  speckle  targets,  speckle  clutter,  and  target 
detection  probability  p^,  *  0.99. 
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be  kept  25-30  dB  below  the  signal  strength  or  the  effective  false  alarm 
probability  will  be  significantly  increased. 

The  average  clutter  return  is  expected  to  be  of  roughly  the  same  strength 
as  the  average  target  return.  Thus,  low  clutter-to-signal  ratios  can  only  be 
achieved  if  the  clutter  and  target  returns  are  well  separated  in  frequency. 
Since  clutter  returns  will  be  distributed  in  frequency  about  zero  Doppler 
shift  and  any  realistic  filter  has  some  small  response  out  to  significant 
distances  from  its  center  frequency,  achieving  a  desirable  false  alarm 
probability  will  impose  a  nonzero  limit  on  the  minimum  detectable  target 
velocity. 

Measurements  of  clutter  frequency  spreads  obtained  at  microwave  frequen¬ 
cies  indicate  that  land  clutter  obeys  a  Gaussian  frequency  distribution12 

2  '1/2  2  2 

p(w)  =  (2tto  )  exp[-(u  / 2o  ]  (38) 

with  standard  deviations  a  between  roughly  0.1  and  1  kmph  depending  on  the 
wind  speed  and  type  of  vegetation.  Since  microwave  radars  look  at  large 
projected  areas  on  the  ground,  one  interpretation  of  these  results  is  that 
the  individual  moving scatterers  (leaves,  branches,  etc.)  have  velocities 
which  are  Gaussianly  distributed.  If  the  angular  resolution  of  the  radar  is 
poor,  many  scatterers  will  be  contained  in  each  pixel  and  each  clutter  return 
will  have  a  broadband  Gaussian  frequency  spread.  If  the  angular  resolution 
is  very  high,  only  a  single  moving  scatterer  will  be  contained  in  each  pixel 
and  a  single-frequency  return  (neglecting  the  linewidth  of  the  incident 
radiation)  will  result.  In  this  limit  Eq.  (38)  gives  the  probability  of 
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obtaining  a  return  at  any  given  Doppler  shift.  The  system  performance  in  these 
two  limits  will  be  markedly  different. 

Let  us  consider  the  broadband  limit  first.  An  indication  of  the  minimum 
detectable  velocity  can  be  be  obtained  by  considering  a  frequency  analyzer 
which  consists  of  a  bandpass  filter  whose  center  frequency  can  be  varied.  We 
assume  a  filter  function  of  the  form 

-1 

a)  -  a)  12 

fCm-OL,)  -  [1  +  (— -)  ]  (39) 

o  tac 

where  u)Q  is  the  center  frequency  and  oic  is  the  3-dB  half -bandwidth  of  the 
filter.  Further  clarification  of  this  system  can  be  found  in  Fig.  13.  Eq. (39) 
approximates  the  response  of  a  6-cavity  LC  filter  of  bandwidth  2wc.  Given  a 
radar  return  with  a  Gaussian  spectrum  centered  at  zero  frequency,  the  signal 
coming  out  of  the  filter  is  proportional  to  the  quantity  I 


OO 

I(V  u>c,  o)  *  £  p(w)  f(u>-<D0)  du  (40) 

We  will  now  make  the  assumption  that  the  target  has  the  same  frequency 
spread  as  the  clutter,  but  is  shifted  in  frequency.  Using  this  assumption, 
which  will  be  shown  shortly  to  be  better  than  it  first  appears,  the  clutter- 
to-signal  ratio  may  be  written  as 


CNR  _  I((V  V  ^ 
SNR  1(0,  u)c,  a) 


(41) 
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Fig.  13.  Functional  relationships  of  MTi  clutter  model.  The  solid  curve 
represents  a  Gaussian  frequency  distribution  with  standard  deviation  a  *  1.0. 
The  dashed  curve  represents  the  response  of  a  6-cavity  LC  filter  of  bandwidth 
2ojc  (*  1.0  a  in  this  exanple)  and  center  frequency  o^(s  5.0  a). 
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If  we  identify  the  Doppler  resolution  with  the  quantity  2u)c  and  the  minimum 
detectable  Doppler  shift  with  the  value  cjq  such  that  CNR/SRN  *  -50  dB  (yield¬ 
ing  an  acceptable  false  alarm  probability) , upon  evaluation  of  the  convolution 
integrals  (Eq.  (40))  we  can  obtain  a  family  of  performance  curves  parametric 
in  a.  In  Fig.  14  the  minimum  detectable  Doppler  shift  is  plotted  as  a  func¬ 
tion  of  the  Doppler  resolution  for  a  =  1.0,  1.25,  1.5,  and  2.0.  These  values 
were  chosen  because  the  actual  frequency  spread  of  clutter  returns  is  deter¬ 
mined  by  the  convolution  of  the  clutter  profile  with  the  laser  frequency  pro¬ 
file  (also  Gaussian) .  Since  at  infrared  frequencies  o(laser)  ^  1  kmph  and 
a(clutter)  ^  0.1  -  1.0  knph,  a  values  between  1  and  2  knph  are  to  be  expected. 

No  units  are  indicated  in  Fig.  14  because  the  amputations  have  been  carried 
out  for  the  case  of  wQ,  u)Q,  and  a  having  the  same  units,  and  are  therefore 
valid  for  ary  consistent  choice  of  units.  Because  o(laser)  dominates  aflut¬ 
ter)  in  many  cases,  target  and  clutter  returns  will  have  roughly  the  same 
frequency  spreads.  It  should  also  be  noted  that  the  oscillations  in  the  curves 
are  real  and  not  artifacts  of  the  conputational  procedure. 

The  effects  of  integrating  returns  from  10  pixels  have  been  estimated  by 
performing  the  same  convolutions  as  before  but  now  identifying  the  minimum 
detectable  velocity  with  the  value  of  such  that  CNR/ SNR  *  -30  dB.  As 
determined  from  Fig.  12,  CNR/SNR  ■  -30  dB  results  in  an  acceptable  effective 
false  alarm  probability  when  10  pixels  are  integrated.  The  results  of  this 
analysis  are  given  in  Fig.  15. 

We  will  now  investigate  the  single-frequency  return  limit.  Referring  to 
Fig.  11  we  note  that  if  a  return  has  CNR/SNR  £  -20  dB,  it  has  a  high  probability 
of  being  detected.  Therefore,  the  effective  false  alarm  probability  (ignoring 
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Fig.  14.  Minimum  detectable  velocity  (Doppler  shift)  versus  Doppler  resolu¬ 
tion  as  a  function  of  the  frequency  width  of  the  returns  (a)  for  the  broadband 
return  Gaussian  frequency  profile  limit  of  clutter  behavior.  The  circles 
represent  numerically  -  calculated  values.  Pp  10"  ,  Ppp  *  0.99. 
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Fig.  15.  Minimum  detectable  velocity  (Doppler  shift)  versus  Doppler  resolu¬ 
tion  as  a  function  of  the  frequency  width  of  the  returns  (a)  for  the  broadband 
return  Gaussian  frequency  profile  limit  of  clutter  behavior.  Returns  from 
10  pixels  have  been  integrated  to  produce  these  results.  The  circles  repre¬ 
sent  numerically  -  calculated  values.  Pi  %  lC^,  p  .  o.99. 


noise  false  alarms)  is  approximately  given  by  the  probability  that  a  single 
frequency  return  occurs  within  the  20  dB  bandwidth  of  the  filter.  This 
probability  is  given  by 


(42) 


(43) 

Perforating  these  integrations  and  identifying  the  minimum  detectable  velocity 

-12 

with  the  valve  of  <i>o,  such  that  the  desired  P'p(%  10  in  this  exanqple)  is 
obtained,  we  obtain  curves  such  as  those  in  Fig.  16. 

The  effects  of  integrating  10  pixels  have  also  been  investigated  in  the 
s ingle -frequency  limit.  Referring  to  Fig.  12  we  note  that  the  detection 
probability  in  the  10-pixel  case  is  roughly  unity  for  CNR/SNR  £  -5  dB.  Thus, 
the  false  alarm  probability  is  approximately  given  by  the  probability  that  a 
single -frequency  return  occurs  within  the  5  dB  bandwidth  of  the  filter. 

This  probability  is  given  by  Eq.  (42)  with  limits  of  integration  of  u>o  +  1.05uic. 
Curves  of  the  minimun  detectable  velocity  as  a  function  of  Doppler  resolution 
for  P'p  %  1012  are  shown  in  Fig.  17. 

The  preceding  clutter  results  are  sraimarized  for  o*  1.0  in  Fig.  18. 
Examination  of  these  curves  indicates  several  features: 

1)  There  is  a  substantial  difference  in  minimum  detectable  velocity 
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Since  we  assune  a  unit  response  when  o)Q  =  0,  we  find 
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Fig.  16.  Minimum  detectable  velocity  (Doppler  shift)  versus  Doppler  resolu¬ 
tion  as  a  function  of  the  width  of  the  probability  distribution  (a)  for  the 
single -frequency  return  Gaussian  probability  distribution  limit  of  clutter 
behavior.  The  circles  represent  numerically  -  calculated  values.  Pi  %  10’  , 
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Fig.  17.  Minimum  detectable  velocity  (Doppler  shift)  versus  Doppler  resolu¬ 
tion  as  a  function  of  the  width  of  the  probability  distribution  (a)  for  the 
single -frequency  return  Gaussian  probability  distribution  limit  of  clutter 
behavior.  Returns  from  10  pixels  have  been  integrated  to  produce  these 
results.  The  circles  represent  numerically  -  calculated  values.  pp  %  10 
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between  the  broadband  and  single -frequency  return  limits.  Since 
infrared  radar  returns  are  expected  to  have  some  of  both  types  of 
characters,  realistic  values  will  lie  somewhere  in  between. 

2)  Pixel  integration  improves  system  performance.  However,  for  real¬ 
istic  numbers  of  pixels  integrated,  this  improvement  is  not  very 
dramatic.  Consequently,  pixel  integration  will  not  be  a  total 
panacea  for  the  minimum  detectable  velocity  problem. 

3)  The  minimum  detectable  velocity  increases  with  increasing  filter 
bandwidth.  Thus,  the  best  performance  is  obtained  with  the  smallest 
filter  bandwidth  compatible  with  the  linewidth  of  the  target  returns. 
However,  the  smaller  the  bandwidth,  the  more  filters  required  to 
cover  the  entire  range  of  expected  target  velocities.  As  a  result, 
the  trade-off  between  performance  and  complexity  must  be  considered. 

4)  For  typical  values  (o  %  2u>c  %  1  knph)  minimum  detectable  velocities 
will  be  of  the  order  of  5-10  knph.  Although  this  is  not  the  1-2  kmph 
mentioned  as  a  design  goal  in  Section  1,  5-10  knph  is  acceptable. 

All  of  the  preceding  results  have  assumed  a  Gaussian  clutter  spectrum. 
However,  a  Gaussian  distribution  has  tails  extending  to  infinity  while  real 
clutter  sources  cannot  have  components  with  infinite  velocities.  In  addition, 
the  experimental  results  quoted  were  obtained  under  conditions  where  large 
numbers  of  independent  scatters  are  contained  within  the  beamwidth.  Under 
such  conditions  the  law  of  large  nunbers  applies  and  the  measured  distribution 
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will  tend  to  a  Gaussian  regardless  of  the  shape  of  the  actual  distribution. 

As  a  result,  the  clutter  spectrum  may  not  and  in  fact  cannot  be  truly 
Gaussian.  Therefore,  Figs.  14-18  should  be  interpreted  as  giving  qualitative 
trends  as  opposed  to  quantitative  results.  For  exanple,  a  sharp  cutoff  of 
the  high-velocity  components  would  tend  to  reduce  the  minimum  detectable 
velocity.  It  is  also  possible  that  the  clutter  distribution  may  be  Gaussian 
out  to  several  standard  deviations.  If  this  is  true,  and  if  the  center 
frequency  of  the  filter  lies  within  the  Gaussian  regime,  then  Figs.  14-18 
will  be  quantitative.  Experimental  investigations  of  the  clutter  spectrum 
at  infrared  frequencies  are  needed  before  the  effects  of  clutter  on  MTI 
performance  can  be  conclusively  evaluated. 

5.  ARCHITECTURE  OF  THE  MTI  ELECTRONICS 

The  basic  corrponents  of  any  airborne  laser  Doppler  MTI  system  are  indi¬ 
cated  schematically  in  Fig.  19.  A  laser  beam  of  known  frequency  content  is 
transmitted  to  the  ground  and  the  radiation  reflected  from  the  target  is  col¬ 
lected  and  detected  by  a  detector  (or  possibly  an  array  of  detectors) .  The 
electrical  signal  from  each  detector  is  then  frequency  analyzed  and  the  tar¬ 
get  return  frequency  data  are  processed  along  with  data  concerning  the  aircraft 
velocity  to  determine  the  presence  or  absence  of  objects  moving  with  respect 
to  the  ground.  Finally,  the  MTI  signals  are  provided  to  a  display  or  to  an 
automatic  cuing  system. 

The  MTI  processing  electronics  must  be  capable  of  performing  several 
operations.  As  shown  in  the  preceding  sections  both  postdetection  integration 
and  multiple  pixel  integration  are  required  for  maximum  signal-to-noise 
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Fig.  19. 


Schematic  diagram  of  an  airborne  laser  Doppler  MTI  system. 
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performance.  In  addition,  if  the  aircraft  lacks  an  adequate  inertial 
platform,  the  MTI  processor  may  be  required  to  generate  its  own  aircraft 
velocity  signal.  Finally,  the  processor  must  incorporate  some  means  of 
using  the  target  velocity  data  to  indicate  moving  targets  to  the  pilot. 

This  last  function  can  be  accomplished  in  two  basically  different 
ways.  In  the  first,  the  strongest  velocity  component  coming  from  each 
pixel -integrated  region  of  the  field-of-view  is  determined  and  used  to 
produce  a  color-coded  image  (different  velocities  would  be  coded  in  dif¬ 
ferent  colors  in  much  the  same  manner  as  is  done  in  range/intensity  color 
mapping  )  which  is  rolled  down  a  display  as  the  aircraft  flies  along. 

The  pilot  would  examine  this  display  and  manually  switch  into  the  imager 
mode  to  investigate  suspected  targets.  In  an  alternate  approach  the 
pixel -integrated,  frequency- analyzed  data  would  be  notch- filtered  (to 
eliminate  clutter  returns)  and  threshold  detected.  Any  signal  exceed¬ 
ing  the  threshold  would  cause  the  infrared  radar  to  automatically  switch 
into  the  imaging  mode  and  start  generating  images  of  the  region  immedi¬ 
ately  surrounding  the  pixels  which  exceeded  the  threshold.  Since  the 
pilot  will  in  general  be  too  busy  flying  the  airplane  to  continually  scan 
an  MTI  display,  we  will  opt  for  the  latter  approach  because  of  its  auto¬ 
matic  cuing  character. 

Let  us  now  consider  the  problem  of  determining  the  aircraft  velocity. 

One  of  the  simplest  approaches  is  found  in  TACCAR  (short  for  time-averaged- 

12 

clutter  coherent  airborne  radar)  systems.  A  schematic  diagram  of  a 
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TACCAR-like  velocity  correction  system  is  shown  in  Fig.  20.  The 
heterodyne  signal  from  a  detector  passes  through  a  mixer  and  a 
portion  of  the  mixer  output  feeds  a  frequency  discriminator.  The 
discriminator's  output  controls  a  frequency  synthesizer  which  provides 
the  second  input  to  the  mixer.  This  servo  loop  attempts  to  maintain 
the  mixer  output  frequency  at  a  constant  value  wQ.  However,  the  tem¬ 
poral  response  of  the  servo  loop  is  purposefully  made  slower  than  the 
characteristic  time  scales  over  which  the  input  frequency  can  change 
(i.e.,  the  dwell  time)  but  fast  compared  to  the  time  scales  over  which 
the  aircraft  velocity  can  significantly  change.  Thus,  the  servo  system 
will  track  out  changes  in  signal  frequency  due  to  aircraft  (and  scanner) 
motions  but  will  not  respond  to  rapid  frequency  changes  associated  with 
the  appearance  of  a  moving  target  in  the  laser  beam.  In  these  systems 
the  output  of  the  mixer  is  typically  passed  through  a  band-reject  filter 
which  eliminates  frequencies  close  to  u)Q.  Signals  from  stationary  objects 
(which  will  occur  at  ojq)  will  be  rejected  while  the  signals  from  moving 
targets  (occurring  at  frequencies  u>  f  u)Q)  will  be  passed.  The  moving 
target  signals  are  then  detected,  integrated  to  reduce  noise  (the 
total  passband  of  the  band- reject  filter  must  be  large  compared  to 
the  dwell  bandwidth  in  order  to  allow  all  non-zero  velocities  of 
interest  to  be  detected),  and  compared  against  a  threshold.  Signals 
exceeding  the  threshold  are  identified  as  moving  targets  and  sent  to 
the  display  or  automatic  cuing  system. 


Fig.  20.  Schematic  diagram  of  a  TACCAR-like  MTI  system  using  a  servo  loop 
to  conpensate  for  aircraft  and  scanner  motions. 


The  TACCAR  approach  is  not  applicable  to  the  IRAR  MTI  problem  for  two 
reasons  both  of  which  are  traceable  to  the  low  signal -to -noise  ratios 
encountered  in  the  close  air  support  scenario.  First,  TACCAR  frequency 
discriminators  operate  by  splitting  an  input  signal  into  two  parts  and  delay¬ 
ing  one  part  with  respect  to  the  other.  The  two  signals  are  then  combined 
in  a  mixer  which  produces  a  dc  voltage  proportional  to  the  relative  phase 
between  the  two  signals  (and  thus  to  the  frequency  if  a  nondispersive  fixed 
delay  is  used).  The  phase  of  a  signal  can  only  be  well-defined  if  the  signal - 
to-noise  ratio  is  high.  At  the  signal -to -noise  ratios  encountered  in  the 
IRAR  problem,  frequency  discriminators  do  not  perform  satisfactorily.  Second, 
band-reject  filter-integrator-threshold  detector  systems  do  not  perform  well 
at  low  signal-to-noise  ratios.  Consider  the  idealized  band- reject  filter  in 
Fig.  21  with  a  rejection  band  of  width  2W^  centered  about  and  passbands 
of  width  on  either  side  of  the  rejected  band  followed  by  an  integrator 
with  integration  time  T.  If  W2  »  Wj  and  we  consider  ideal  (i.e.,  square) 
filter  and  integrator  responses,  the  signal-to-noise  ratio  of  the  system  is 
given  by  ** 


where  n  is  the  quantum  efficiency,  is  the  incident  power,  and  hv  is  the 
photon  energy.  Noting  that  the  heterodyne  detection  signal-to-noise  ratio 
of  a  system  with  bandwidth  1/T  is 


SN^hT  =  nP^T/hv 


(45) 


Eq.  (44)  may  be  rewritten  as 


SNR  % 


(SNRhT)2 

2W2T 


If  SNR^t/W2T  »  1,  Eq.  (46)  reduces  to 


(46) 


SNR  £  SNRhT/2 


(47) 


while  if  SNRj^/l^T  «  1,  we  find 
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2W2T 


(48) 


Since  typical  values  of  W2T  are  in  the  range  100-1000  and  typical  values  of 
SNR^T  will  be  less  than  100,  Eq.  (48)  indicates  that  the  effective  signal - 
to-noise  ratio  will  be  much  reduced  from  its  theoretical  maximum.  Thus, 
TACCAR-like  systems  are  not  desirable. 

An  MTI  processor  which  does  not  suffer  these  drawbacks  and  incorporates 
all  of  the  necessary  functions  is  shown  in  Fig.  22.  The  output  from  each  of 
Np  detectors  is  frequency  analyzed  into  a  number  (N^)  of  Doppler  channels . 
Several  potential  methods  of  performing  the  frequency  analysis  are  discussed 
in  Section  6.  Each  Doppler  channel  is  individually  envelope  detected  and 
integrated.  The  integration  time  is  chosen  to  correspond  to  several  dwell 
times  (e.g.,  ^500  ysec  if  tD  100  ysec) .  This  accomplishes  both  the  post¬ 
detection  integration  and  one  dimension  of  the  pixel  integration.  Next  all 
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of  the  integrated  outputs  are  now  processed  by  a  multiple  pixel  integrator. 
This  unit  groins  the  Doppler  channels  according  to  frequency  (refer  to  Fig.  23) 
and  for  each  frequency  it  sums  the  signals  from  various  groups  of  detectors. 
This  second  dimension  of  pixel  integration  results  in  M(<  Np)  target  eleva¬ 
tion  channels . 

Each  of  the  outputs  from  the  multiple -pixel  integrator  is  split  into  two 
parts.  In  one  branch  the  outputs  from  the  elevation  channels  of  each  fre¬ 
quency  are  summed  and  tenporally  integrated  with  integration  times  being  many 
times  the  dwell  time  (e.g.,  a.  5  msec).  The  integrated  outputs  (one  for  each 
Doppler  channel)  are  sent  to  a  peak  detector  which  determines  which  channel 
has  the  maximum  signal.  A  disable  signal  is  sent  to  the  disable  switches  in 
the  corresponding  Doppler  channel  in  the  second  branch.  The  characteristic 
response  time  of  the  signal  at  the  input  to  the  peak  detector  in  the  first 
branch  is  much  slower  than  the  response  time  of  the  signal  at  the  output  of 
the  multiple  pixel  integrator.  Thus,  although  the  disable  circuit  will  com¬ 
pensate  for  slow  temporal  variations  in  the  signal  (such  as  those  caused  by 
scanner  motion  and  aircraft  velocity  changes) ,  it  will  not  compensate  for 
fast  variations  (such  as  those  caused  by  the  sudden  appearance  of  a  small 
moving  target).  It  should  be  noted  that  because  of  clutter  effects,  disable 
signals  may  also  have  to  be  sent  to  switches  in  several  Doppler  channels  on 
either  side  of  the  maximum  to  reduce  false  alarms.  This,  of  course,  deter¬ 
mines  the  minimum  detectable  velocity. 

Next,  signals  passed  by  the  disable  switches  (and  thus  which  correspond 
to  non- zero  velocities)  are  threshold  detected.  The  threshold  detector  sig¬ 
nals  are  sent  to  a  peak  detector  which  determines  if  a  moving  target  is  pres¬ 
ent,  which  elevation  channel  has  the  maximun  signal,  and  what  is  the  target 
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OUTPUTS  (M  Elevation  Channels,  Nc  Doppler  Channel*  Per  Elevation  Channel) 


Fig.  23.  Flow  diagram  of  the  multiple  pixel  integrator.  The  upper  portion 
of  the  device  performs  a  sorting  of  the  outputs  according  to  frequency  while 
the  lower  portion  suns  (for  any  single  frequency)  the  outputs  from  different 
gimps  of  detectors.  For  simplicity  only  5  detectors  and  3  Doppler  channels 
per  detector  are  depicted.  A  similar  diagram  could  be  drawn  for  any  desired 
system. 


velocity.  An  enable  signal  from  the  peak  detector  activates  a  switch  which 
converts  the  infrared  radar  from  the  MTI  mode  to  the  imager  mode.  The 
elevation  signal  is  coupled  with  an  azimuth  signal  from  an  encoder  on  the 
pointer-scanner  mirror  to  direct  the  system  to  point  the  imager  in  the  proper 
direction. 

The  MTI  processor  of  Fig.  22  is  clearly  not  unique  and  other  candidate 
systems  can  be  readily  devised.  The  final  design  will  depend  in  part  on  the 
technology  available.  For  example,  the  system  described  above  is  quite  com¬ 
plicated  (as  are  other  candidate  systems)  and  has  the  potential  for  becoming 
quite  large.  However,  many  of  the  necessary  operations  can  be  performed 
either  in  a  digital  or  discrete  analog  fashion.  It  is  envisioned  that 
charge -coup led  device,  microprocessor,  and  very  large-scale  integration 
(VLSI)  technologies  may  be  potentially  useful  in  reducing  the  size  of  the 
processing  electronics.  The  nature  (e.g.,  serial  or  parallel  outputs  and 
inputs)  of  the  devices  employed  will  clearly  affect  the  details  of  the  design. 
Nevertheless,  any  proposed  system  must  incorporate  most  of  the  features  of 
Fig.  22. 

6.  FREQUENCY  ANALYSIS  TECHNIQUES 

There  are  three  major  candidate  techniques  for  performing  the  frequency 
analysis  function:  bandpass  filter  banks,  acoustooptic  deflection,  and  chirp 
z  transform  Fourier  analysis.  The  first  technique  involves  nothing  more  than 
dividing  the  amplified  signal  from  the  heterodyne  detector  into  m  ♦  1  parts 
and  passing  each  through  a  different  bandpass  filter  (Fig.  24).  The  filters 
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Fig.  24.  Schematic  diagram  of  a  frequency  analyzer  employing  a  filter  band 


are  chosen  to  form  a  picket  fence  with  partially  overlapping  passbands.  Plac¬ 
ing  the  overlap  at  the  -3  dB  response  points  of  the  filter  gives  a  good  com¬ 
promise  between  channel  ambiguity  (having  a  monochromatic  signal  produce 
outputs  from  several  filters  simultaneously)  and  signal  miss  probability  (not 
detecting  a  weak  monochromatic  signal  which  falls  midway  between  the  center 
frequencies  of  two  adjacent  filters).  The  picket  fence  response  of  an  m  +  1 
filter  bank  with  -3dB  overlap  is  shown  in  Fig.  25.  The  output  from  each  fil¬ 
ter  is  then  envelope  detected  and  supplied  to  the  Doppler  processor.  The  net 
result  of  using  the  filter  bank  is  to  generate  a  discrete  Fourier  transform 
of  the  heterodyne  detector  output.  The  filter  bank  of  Figs.  24  and  25  would 
produce  m  +  1  frequency  components  with  resolution  Aw  covering  an  input 
frequency  range  from  wQ  -  1/2  Aw  to  wQ  +  (m  +  1/2)  Aw.  In  an  MTI  system, 

Aw  is  determined  by  the  required  velocity  resolution  and  the  input  bandwidth 
(and  hence  the  nunber  of  channels  required) 

Bj  *  (m  +  1)  Aw  (49) 

is  determined  by  range  of  Doppler  frequencies  which  may  be  encountered. 
Frequency  analysis  using  bandpass  filter  banks  has  been  a  proven  technique 
for  many  years.  However,  the  microwave  components  are  rather  bulky  (a  100 
channel  filter  bank  might  occupy  3000  -  5000  cm3  and  weigh  5-10  kg)  and 
therefore,  this  technique  may  not  be  practical  for  an  airborne  multichannel, 
multidetector  Doppler  analyzer.  Acoustooptic  and  chirp  z  transform  devices, 
on  the  other  hand,  offer  significant  potential  savings  in  size  and  weight 
when  high  resolution,  wide  bandwidth  devices  are  required. 
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RELATIVE  RESPONSE  (48) 
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Fig.  25.  Picket  fence  response  of  a  filter  bank  of  m  +  1  filters  which  over¬ 
lap  at  the  -3  dB  response  points. 
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Acoustooptic  frequency  analyzers  work  on  the  principle  that  the  angular 
deflection  that  an  optical  beam  suffers  when  it  interacts  with  an  acoustic 
wave  is  a  function  of  the  frequency  of  the  acoustic  wave.  As  indicated  in 
Fig.  26,  in  an  acoustooptic  frequency  analyzer  the  signal  from  the  heterodyne 
detector  is  anplified  and  transformed  into  an  acoustic  wave  in  a  crystal 
(such  as  LiNbOj)  by  a  transducer.  A  fraction  of  the  intensity  of  a  visible 
laser  beam  incident  on  the  crystal  is  deflected  through  various  angles  related 
to  the  frequency  components  of  the  acoustic  wave  and  focused  onto  a  detector 
array.  Since  the  intensity  of  the  laser  radiation  deflected  through  any  given 
angle  is  proportional  to  the  amplitude  of  the  corresponding  acoustic  frequency 
component,  the  magnitudes  of  the  detector  outputs  constitute  a  discrete 
Fourier  transform  of  the  heterodyne  detector  signal.  Fig.  26  illustrates 
this  behavior  for  an  input  signal  containing  two  discrete  frequencies  fj  and 

f2* 

The  theory  of  acoustooptic  frequency  analysis  has  been  thoroughly  reviewed 
14  15 

in  the  literature.  ’  However,  for  completeness  we  will  sunmarize  the  key 
results.  Consider  an  acoustic  grating  of  wavelength  A  formed  in  a  crystal 
by  an  acoustic  wave  of  frequency  f  with 

A  =  V/f  (50) 

where  V  is  the  sound  velocity.  A  light  beam  of  free  space  wavelength  XQ 
incident  tpon  this  grating  at  the  Bragg  angle  0£  given  by 

sin  e£  »  K/2k  (51) 
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COMPONENT 

OUTPUTS 


Fig.  26.  Schematic  diagram  of  an  acoustooptic  frequency  analyzer.  The 
acoustic  wave  patterns  correspond  to  an  input  containing  two  frequencies 
f^  and  $2' 
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where  K  =  2ir/A,  k  =  2nn/  XQ,  and  n  is  the  refractive  index,  will  be  diffracted. 
Conservation  of  mcmentun  (Fig.  27a)  requires  the  diffracted  beam  angle  0d  to 
equal  the  Bragg  angle  (i.e.,  the  angle  of  incidence  0!) 

ed  “  8B  =  6i 

In  this  treatment,  primed  angles  are  measured  in  the  crystal  while  unprimed 
angles  are  measured  external  to  the  crystal.  The  two  sets  of  angles  are 
related  by  Snell's  law 

n  sin  0*  =  sin  0  (52) 


If  for  a  fixed  angle  of  incidence  0|,  the  acoustic  wavelength  is  close  but 
not  equal  to  the  wavelength  at  which  the  Bragg  condition  is  satisfied,  a 
diffracted  beam  will  still  result  even  though  momentum  cannot  be  strictly 
conserved.  The  diffracted  beam  will  occur  at  the  angle  0^  which  minimizes 
the  momentum  mismatch.  Analysis  of  Fig.  27b  yields 

tan  0^  *  (K  -  ksin  0!)  /  kcos  0j  (53) 


The  momentum  mismatch,  of  course,  will  result  in  a  reduced  diffracted  intensity. 
The  diffracted  intensity  normalized  to  the  incident  intensity  may  be  expressed 
as 

Id/I  %  [sin2  q1/2]  [sinc2(-L-  (F2  -  F))]  (54) 

2L0 

n  sine2  [~  (F2  -  F)]  (55) 

2Lo 
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Fig.  27.  Momentum  conservation  in  the  acoustooptic  interaction,  (a)  the  Bragg 
condition  and  (b)  nomentun  mismatch  at  non-optimum  angles  of  incidence. 


where  F  »  f/fQ  is  the  acoustic  frequency  normalized  to  the  frequency  at 
which  Bragg  matching  occurs,  L  is  the  width  of  the  acoustic  beam. 


L0  =  nA0  cos6!/X0 


and  Aq  =  V/fQ.  The  diffraction  efficiency  n  is  given  by 


-4(*¥)- 

2\  \  pV*  /  H 


cos  0!^ 


where  p  is  the  crystal  elastooptic  coefficient,  p  is  the  material  density, 
and  Pg  is  the  acoustic  power  in  a  beam  of  cross  sectional  area  L  •  H.  From 
Eq.  (55)  the  3  dB  bandwidth  of  the  acoustooptic  frequency  analyzer  is  found 


to  be 


Af3dB  ^  1-8  foLo^L 


The  frequency  resolution  <5f  of  an  acoustooptic  frequency  analyzer  is  related 
to  the  nunber  of  acoustic  wavelengths  N  encompassed  by  the  laser  beam.  If 
D  is  the  projected  width  of  the  laser  beam,  the  resolving  power  is  given  by 

R  -  X/6X  =  f/6f  =  N  =  D/A  *  Df/V  (59) 


6f  -  V/D  -  1/x 


where  t  is  the  acoustic  transit  time  across  the  laser  beam.  The  total  number 
of  resolvable  frequency  elements  N  contained  in  the  bandwidth  of  the  device 
is  given  by  the  time -bandwidth  product 

N  *  ^3  dB  /  =  Af3  dB  T  *-61^ 

Both  bulk-wave  and  surface-wave  acousto-optic  frequency  analyzers  have 
been  fabricated  with  time -bandwidth  products  in  excess  of  1000. 14  16  The 
latter,  using  integrated  optics  (guided  optical  waves,  surface  acoustic  pro¬ 
cessing),  offers  potential  for  significant  miniaturization.  Unfortunately, 
the  integrated  optical  devices  probably  cannot  achieve  the  desired  frequency 
resolution.  This  results  because  typical  acoustic  velocities  are  of  the 

c 

order  of  a  few  times  10  cm/sec  implying  laser  beam  widths  of  the  order  of  a 
few  cm  to  achieve  100  kHz  resolution.  Bulk  wave  devices  can  be  made  this 
large  but  at  the  present  time  integrated  optical  structures  cannot. 

17 

Chirp  z  transform  devices  work  on  the  following  mathematical  principles, 

:  o_r* 

Consider  an  input  signal  S^(t)e  passing  through  a  chirp  transversal 
filter  with  inpulse  response  h(t)  *  ei2Trfot  giat  .  Suppressing  the  carrier 
frequency  dependence,  the  output  of  the  transversal  filter  is  given  by  the 
convolution  of  with  h 

OO  2 

So(T)-/Si(t)  eia(t_T)  dt  (62) 


73 


Expanding  the  exponent  yields 


S0(t)  >  e1”2 /’siCt)eic,t2  e'l2“Tt  dt  (63) 

-o© 

±  2 

which  except  for  the  e10tT  weighting  factor  is  the  Fourier  transform  (with 

iat-2 

transform  variable  u  =  2ax)  of  the  function  S^(t)e  .  If  we  define  a  sig- 

_ 

nal  S’ (t) ,  premultiplied  by  a  chirp  e  such  that 


S^t) 


S!  (t)  e 


(64) 


and  an  output  signal  which  is  the  transversal  filter  output  multiplied  by  a 
chirp 

(x)  *  Sq(t)  e'120"2  (65) 

we  find 

oo 

s-  (t)  S!  (t)  e”2aTt  dt  (66) 

-OO 

That  is ,  (x)  is  the  Fourier  transform  of  Sj  (t) .  If  the  chirps  are  of 

finite  duration  and  repetitive,  a  finite  Fourier  transform  operation  is  per¬ 
formed. 

The  chirp  z  Fourier  transformation  can  be  implemented  by  a  system  shown 

17 

schematically  in  Fig.  28.  Impulse  inputs  to  two  chirp  transversal  filters 

_iat2  j  2 

generate  the  chirps  e  and  e  which  multiply  the  input  and  output 
signals,  respectively.  A  third  chirp  transversal  filter  performs  the  required 
convolution. 
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Fig.  28. 


Schematic  diagram  of  a  chirp  z  transformation  frequency  analyzer. 
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A  chirp  transversal  filter  can  be  considered  as  an  ideal  linear  dis¬ 
persive  delay  line.  In  such  a  delay  line  the  delay  encountered  by  a  mono¬ 
chromatic  wave  is  proportional  to  its  frequency.  Clearly,  the  impulse 
response  of  such  a  linear  dispersive  delay  line  is  a  linear  frequency  chirp. 
The  direction  (sign)  and  rate  of  the  chirp  is  proportional  to  the  frequency 
dispersion  while  the  maximum  amplitude  of  the  chirp  is  given  by  the  bandwidth 
of  the  delay  line. 

Surface  acoustic  wave  (SAW)  devices  can  be  used  to  realize  a  linear  dis¬ 
persive  delay  line  at  microwave  frequencies.  Consider  the  reflective  array 
compressor  (RAC)  shown  in  Fig.  29.  Surface  acoustic  waves  are  generated 
by  a  transducer  and  propagate  through  an  array  of  parallel  grooves  etched  in 
the  surface  of  the  crystal.  The  grooves  are  oriented  at  an  angle  0  to  the 
direction  of  propagation  of  the  incident  surface  acoustic  waves.  When  the 
waves  encounter  a  sequence  of  grooves  whose  spacing  G  is  such  that  construc¬ 
tive  interference  can  occur,  the  acoustic  waves  will  be  reflected  at  right 
angles.  The  depth  of  the  grooves  controls  the  reflection  coefficient.  The 
reflected  waves  encounter  a  second  mirror  image  array,  are  reflected  again, 
and  propagate  to  an  output  transducer.  Since  the  acoustic  velocity  is 
independent  of  frequency,  the  delay  encountered  by  any  acoustic  frequency 
is  proportional  to  the  total  path  length  travelled.  Consequently,  by  control¬ 
ling  the  sequence  of  groove  spacings  different  delays  can  be  obtained  for 
different  frequencies.  Minor  groove  spacing  (phase)  errors  can  be  conpensated 
by  depositing  a  film  of  varying  width  between  the  two  arrays  (as  shown  in 

Fig.  29). 
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In  many  materials,  the  surface  wave  velocities  differ  along  the  incident 
and  transverse  propagation  directions.  Therefore,  the  groove  angle  0  need  not 
be  45°.  Consider  Fig.  30.  If  v^  and  v 2  represent  the  surface  wave  velocities 
in  the  incident  and  transverse  directions,  then  the  requirement  that  reflec¬ 
tion  generates  a  flat  wavefront  at  right  angles  requires  the  path  lengths 
ABD  and  Ac  to  differ  by  a  wavelength, 

JS  +  EC  =  AC  +  Aj^  (67) 

where  ^  2  s  vi  2^'  T^ie  same  Pat^  length  condition  applied  to  the  time 
reversed  wavefront  requires 


OB'  +  SA  *  DC  +  X2  (68) 

Since  A1  =  G/cos  9,  51  =  G/sin  9,  A5  =  15,  and  55  =  1A,  we  have 

G  *  cos  9  =  X2  sin  6  (69) 

from  which  we  determine 

tan  0  *  Xj/X2  *  Vj/v2  (7°) 


A  linear  chirp  will  be  generated  when  the  groove  spacing  is  linearly 
proportional  to  the  distance  from  the  transducer.  Although  it  will  not  be 
demonstrated  here,  the  proper  groove  positions  are  given  by 


nX^ 

2 


(71) 
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where  xn  is  the  position  of  the  n^  groove.  The  inpulse  response  of  this 
filter  is 


i2irf  t  (1  +  Bt/2f  T) 
h(t)  =  e  0  0 


T 

T 


<  t  < 


T 

I 


(72) 


where  f  =  Vj/ and  the  instantaneous  frequency  f  is 

f  =  fQ  +  (73) 

and  varies  linearly  from  f  -  B/2  to  fQ  +  B/2. 

Both  the  acoustooptic  and  the  chirp  z  transform  techniques  are  potentially 
capable  of  accomplishing  the  frequency  analysis  required  by  the  IRAR  MTI  sub¬ 
system  in  packages  of  reasonable  size.  However,  both  technologies  are  still 
developing  and  devices  capable  of  meeting  our  specifications  while  operating 
in  an  aircraft  environment  have  not  yet  been  demonstrated.  Further  work  on 
these  devices  is  clearly  required  before  they  can  be  incorporated  into  an 
airborne  system. 

7.  CONCLUSIONS 

The  work  presented  in  this  report  appertains  to  a  very  specific  system. 
Therefore,  minor  changes  in  initial  design  parameters  can  have  a  significant 
impact  on  the  numerical  conclusions  arrived  at  here.  However,  certain  general 
conclusions  can  be  drawn  which  pertain  to  airborne  infrared  radar  MTI  systems 
as  a  whole. 

First,  weather  penetration  effects,  compactness ,  efficiency,  high  power, 
and  high  frequency  stability  favor  using  a  cw  laser  transmitter  and  hetero- 
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dyne  detection  in  a  Doppler  MTI  system. 

Second,  the  ability  to  completely  cover  a  wide  field-of-view  in  a  target 
acquisition  mode  will  probably  require  the  use  of  large  detector  arrays.  In 
a  line  scan  search  mode  the  optimum  signal-to-noise  ratio  is  achieved  by  a 
linear  detector  array  with  just  enough  detectors  to  adequately  cover  the 
field-of-view. 

Third,  pixel  integration  will  be  required  in  many  situations  to  produce 
desired  detection  and  false  alarm  probabilities.  In  any  situation  pixel  inte¬ 
gration  will  improve  the  weather  penetration  capability  of  a  system. 

Fourth,  ground  clutter  will  probably  limit  obtainable  minimum  detectable 
target  velocities  to  values  of  the  order  of  5  kmph  or  larger  even  though  the 
velocity  resolution  may  be  of  the  order  of  1  kmph. 

Fifth,  the  size  of  the  MTI  electronics  can  impose  a  serious  constraint 
on  system  utility.  However,  careful  architecture  of  the  processing  electronics 
as  well  as  judicious  use  of  developing  acoustooptic,  charge-coupled  device, 
very  large-scale  integration,  and  surface  acoustic  wave  technologies  may 
ameliorate  this  situation. 

Lastly,  an  MTI  system  can  share  most  of  the  components  of  an  imaging 
infrared  radar  system.  This  should  enable  multifunction  infrared  radar  systems 
to  be  built  using  a  single  optical  train. 
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